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ABSTRACT 

We report the detection of eclipses in GJ 3236, a bright (/ = 11.6) very low mass binary system 
with an orbital period of 0.77 days. Analysis of light- and radial velocity curves of the system yielded 
component masses of 0.38 ± 0.02 and 0.28 ± 0.02 M©. The central values for the stellar radii are 
larger than the theoretical models predict for these masses, in agreement with the results for existing 
eclipsing binaries, although the present 5% observational uncertainties limit the significance of the 
larger radii to approximately la. Degeneracies in the light curve models resulting from the unknown 
configuration of surface spots on the components of GJ 3236 currently dominate the uncertainties in 
the radii, and could be reduced by obtaining precise, multi-band photometry covering the full orbital 
period. The system appears to be tidally synchronized and shows signs of high activity levels as 
expected for such a short orbital period, evidenced by strong Ha emission lines in the spectra of both 
components. These observations probe an important region of mass-radius parameter space around 
the predicted transition to fuUy-convective stellar interiors, where there are a limited number of precise 
measurements available in the literature. 

Subject headings: binaries: eclipsing - stars: low-mass, brown dwarfs - stars: individual (GJ 3236) 



1. INTRODUCTION 

Detached, double-lined eclipsing binaries provide a 
largely model-independent means to precisely and accu- 
rately measure fundamental stellar properties, particu- 
larly masses and radii. In the best-observed systems the 
precision of these can be at the < 1 per cent level, and 
thus pla ce stringent con straints on stellar evolution mod- 
els (e.g. lAndersen|[l991[ ). 

Despite this, as far as we are aware, there are at 
present only four known systems with one or more com- 
ponents between 0.4 and the hydrogen burning limit 
on the main sequence (i .e. old, field st ars): CM Dra 
jEggen & Sandaee 'l967l iLacvl Il977t 'Metcalf e etal 
1996,: Morales et al. 20Qi), CU Cnc B f.Delfosse et al . 
19991 : iRibasI I2003L LP 133- 373 (iVaccaro et al.l I2007Il 



and SDSS J031 824-010018 (iBlake et al.l I2008D. Al- 
though JW 380 (llrwin etldl l2007aD . NSVS 02502726B 
(iCakirh. Ibanoglu fc GunTOJ l2009l. and the NGC 1647 
svstem of iHebb et al. I (I2007D also satisfy this mass crite- 
rion, these objects are still contracting on the pre-main- 
sequence. 

Of the known systems, the only two with parameters 
determined to better than 2 per cent are CM Dra and 
CU Cnc. For SDSS J031824-010018, the knowledge of 
the parameters is limited largely by radial velocity er- 
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rors, since this system is extremely faint (SDSS r = 19.3) 
and has an extremely short orbital period (0.41 days), 
meaning long integrations cannot be used to obtain bet- 
ter signal-to-noisc. Hence there is little possibility for 
substantial improvement in the parameters in the near 
future. It is clear that, in order to better-constrain the 
stellar mass-radius relation on the main sequence, more 
bright, low-mass eclipsing binaries are needed to yield 
extremely precise masses and radii. 

Furthermore, observations of these systems have in- 
dicated significant discrepancies with the stellar mod- 
els. This is particularly the case below 0.4 Mq, and the 
components of CM Dra (the lowest-mass system with 
< 1% observational errors) have radii 10 — 15 per cent 
larger than the theoretical predictions from state-of-the- 
art stellar evolution models. It has been suggested (e.g. 
IChabrier et al.ll2007h that the reason for this discrepancy 
may be that close binaries are not actually representa- 
tive of single-stars at this level of precision. The effect 
of the close companion and tidal locking is likely to sig- 
nificantly increase activity levels in close binaries, and 
these authors suggest that it could be this effect that is 
r esponsible for the infla ted radii of CM Dra. 

iLopez-MoralesI (|2007t ) examined the available sample of 
single-star and eclipsing binary measurements to search 
for correlations of radius with activity levels and metal- 
licity, finding that such a correlation of activity with ra- 
dius does appear to exist for members of close binary 
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stars. This result was based on a small number of mea- 
surements for late M-dwarfs, and the conclusions would 
be strengthened by the availability of additional precise 
estimates of the parameters for binaries with a range of 
activity levels and orbital periods, to explore the avail- 
able parameter space. 

We present the discovery of a new bright, low-mass 
eclipsing binary system, GJ 3236. The components show 
high activity levels, as evidenced by Ha emission lines in 
the spectra of both stars, and show rotational modula- 
tions in the light curves that are synchronized with the 
binary orbital period of 0.77 days. This system thus has 
the potential to yield an additional precise test of stel- 
lar evolution models for high activity and short orbital 
periods. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. MEarth photometry 

Eclipses in GJ 3236 were detected during 2008 January 
in the first two weeks of routine operations of the newly- 
commissioned MEarth observatory, a system designed 
primarily to search for transiting super-Earth exoplanets 
orbiting around the nearest 2000 mid to late M-dwarfs 
in the northern hemisp here (jNutzman fc Charbonneaul 
I2008t llrwin et al.l 120081 ). Exposure times on each field 
observed by MEarth are tailored to achieve sensitivity to 
a particular planet size for the assumed stellar parame- 
ters of the target star, and were 82 s for GJ 3236. 

After the initial detection of the eclipses, we switched 
to a follow-up mode, observing at the highest possible ca- 
dence (i.e. continuously, resulting in a cadence of approx- 
imately 2 minutes including overheads) during eclipse 
and for an additional 1 hour window surrounding the 
eclipse, and at the normal cadence for the MEarth sur- 
vey of approximately 20 minutes for the remainder of the 
time. In this way, we sample as well as possible the rapid 
flux decrement during eclipse, and the out of eclipse por- 
tions of the light curve. We note that sampling the latter 
is of crucial importance for obtaining accurate parame- 
ters from light curve analyses of binary systems showing 
out of eclipse modulations, such as the present one. 

We obtained data both in the initial season after the 
eclipses were first detected (2008 January to May, inclu- 
sive), and following the summer monsoon (2008 October 
onwards). During this time, there were a number of soft- 
ware improvements (including the ability to update the 
telescope pointing after taking each science image, result- 
ing in substantially improved light curves due to mini- 
mizing the drift across the detector, and hence the effect 
of flat flelding errors, fringing, etc.) and our filter sys- 
tem was changed from Cousins / to a custom long-pass 
fllter with transmission from 715 nm, and limited at the 
long wavelength e nd by the tail of the CCD quantu m ef- 
ficiency curve fsee lNutzman fc Charbonneaul l2008l ). We 
have elected to use only the latter data-set for modeling, 
despite this having the disadvantage of being observed 
in a non-standard bandpass, because it is approximately 
contemporaneous with our 1^-band follow-up photome- 
try. We see evidence for evolution of the amplitude and 
phase of the out of eclipse modulation between the two 
observing seasons, so this simplifies the analysis by allow- 
ing the use of a single set of spot parameters to describe 
both bandpasses. 

A total of 1540 observations were obtained between 



UT 2008 October 2 and UT 2008 December 11, includ- 
ing 6 primary and 3 secondary eclipses sampled at high 
cadence. 

Data were reduced using the standard MEarth reduc- 
tion pipeline, which is at prese nt identical t o the M oni- 
tor project pipeline described in ilrwin et al.l (|2007bh . We 
used a smaller aperture radius of 5 pixels (3'.'8) than the 
usual value of 10 pixels adopted for bright stars, due 
to the presence of a nearby star approximately 7" from 
GJ 3236 at the epoch of our observations, and 4 mag 
fainter in the MEarth bandpass. This star does not share 
a common proper motion with GJ 3236 and is therefore 
not physically associated. Reducing the aperture size in 
this way reduces the flux contributed by this star in the 
light curve of GJ 3236 to a negligible level (we estimate 
a 0.1% error in the measured fluxes, and thus also in the 
eclipse depths). 

The high proper motion of the GJ 3236 system allows 
us to constrain the contribution of any additional back- 
ground stars in the photometric aperture that are not 
co-moving with GJ 3236 itself, using previous epochs of 
imaging. We show in Figure [T] a series of three images 
centered on the position of the photometric aperture used 
for GJ 3236 in the MEarth images, demonstrating that 
there are no such objects. If any additional third light 
is indeed present, it must therefore be co- moving with 
GJ 3236, meaning it is highly likely to be physically as- 
sociated, given the extremely small probability of such a 
chance alignment. 

The photometry was calibrated using observations of 
standard star fields from Landolt (|1992f ) taken automat- 
ically each night by the MEarth observation scheduling 
software. By fitting all the standard star observations 
from several photometric nights, we derived the following 
color equation to convert between the MEarth bandpass 
and the standard Johnson-Cousins system: 

MEarth = /c - 0.18 (^j - Ic) (1) 

The full MEarth light curve is reproduced in Table [T] 

2.2. FLWO 1.2 m V -band follow-up photometry 

Observations centered around the primary eclipse of 
UT 2008 September 25 and the secondary echpse of UT 
2008 September 27 were obtained using the KeplerCam 
instrument on the FLWO 1.2 m telescope. We used the 
standard binning 2x2 readout mode, since the plate scale 
of 0'.'34 per unbinned pixel significantly oversamples the 
typical seeing at FLWO. The resulting plate scale was 
0'.'67 per summed pixel. We used the V filter and an 
exposure time of 120 s. We also used observations of 
standard star fields from Landolt (1992) taken on UT 
2008 September 28 to calibrate the KeplerCam photom- 
etry onto the standard Johnson-Cousins system. A total 
of 245 observations were taken, where possible starting 
2 hours before first contact and finishing 2 hours after last 
contact to sample the out-of-eclipse portions of the light 
curve and thus allow the eclipse depths to be properly 
measured. The end of the primary eclipse was missed 
due to twilight. 

These photometric data were reduced using the same 
pipeline as described in ^2.1\ We used an aperture radius 
of 5 binned pixels, corresponding to 3'.'4 on-sky (recall- 
ing from H2.1\ that the size of this aperture is important 
due to the nearby star). The difference in magnitude 
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TABLE 1 
MEarth light curve of GJ 3236. 



HJD 


MEarth 


Error^ 


Am'' 


FWHM (pix)': 


Airmass 


X (pix)'' 


y (pix)'* 


2454741.691783 


11.1329 


0.0036 


-0.278 


7.24 


1.94119 


1029.66 


1048.57 


2454741.704376 


11.1151 


0.0035 


-0.337 


7.18 


1.85590 


1030.44 


1049.71 


2454741.716182 


11.1159 


0.0035 


-0.223 


6.78 


1.78328 


1030.33 


1048.66 


2454741.727410 


11.1127 


0.0035 


-0.209 


6.81 


1.72032 


1030.19 


1048.31 


2454741.738267 


11.1180 


0.0035 


-0.226 


6.95 


1.66467 


1030.41 


1048.00 



Note. — Tablc[l]is published in its entirety in the electronic edition of the Astrophysical Journal. 
A portion is shown here for guidance regarding its form and content. 

Estimated using a standard CCD noise model, including contributions from Poisson noise in the 
stellar counts, sky noise, readout noise and errors in the sky background estimation. Correction 
to the frame magnitude zero-point applied in the differential photometry procedure. More negative 
numbers indicate greater losses. Median FWHM of the stellar images on the frame. The plate scale 
was 0'.'757/pix.<' x and y pixel coordinates of GJ 3236 on the CCD image, derived using a standard 
intensity-weighted moments analysis. 



between GJ 3236 and the fainter star is only 2.3 mag in 
V^-band, which leads to approximately a 0.6% error in 
the measured fluxes of GJ 3236, and thus a similar error 
in the measured eclipse depths. Since the point spread 
functions on our images are not very well-behaved it will 
be challenging to reduce this using PSF-fitting photom- 
etry, and it is still substantially smaller than the other 
sources of error in the present light curve models, so we 
have not pursued this issue further. 
The full F-band light curve is reproduced in Tabled 

2.3. PAIRITEL J -band follow-up photometry 

Observations in the near-infrared J-band were ob- 
tained using the robotic Peters automated infrared imag- 
ing telescope (PAIRITEL) from UT 2008 February 17 to 
UT 2008 March 3 (inclusive). We scheduled observa- 
tions around the times of primary and secondary eclipse, 
obtaining 457 data points spread over 10 nights. Ex- 
posure times were 7.8 s. The observations were auto- 
matically scheduled, collected, and r educed by the fully 
robotic PAIRITE L observing system (jBloom et al.ll2006l : 
iBlake et"aI]|2005D . We then produced differential pho- 
tometry using a set of comparison stars chosen from the 
2MASS catalog. We estimate a photometric error of ap- 
proximately 2% per data point from the scatter of the 
out-of-eclipse measurements. The full J-band light curve 
is reproduced in Table [3] 

2.4. Spectroscopy 

Spectroscopic observations were obtained using the 
TRES fiber- fed echelle spectrograph on the FLWO 1.5 m 
Tillinghast reflector. We used the medium fiber (2'.'3 pro- 
jected diameter) throughout, yielding a resolving power 
of i?~ 48 000. 

Table H] summarizes these measurements, including the 
radial velocities of both components of the binary derived 
therefrom. The TRES instrument is extremely stable, so 
we acquired ThAr wavelength calibration exposures be- 
fore or after the target exposures, rather than using the 
simultaneous calibration fiber, which can lead to con- 
tamination of the target spectrum in the red-most orders 
from the very strong Ar lines in the ThAr lamp spectrum. 
A second fiber was placed on sky, but this was not used 
due to a problem causing extremely poor throughput in 
the TRES sky fiber, which has since been resolved. 



Spectr a were redu ced using standard procedures in 
IRAF^ ()Todvl [T993f) from the echelle package to ex- 
tract the spectra from the target fiber and sky fiber si- 
multaneously. The TRES CCD is read out using two am- 
plifiers, which necessitates combining them before trac- 
ing and extraction. This was done using the mscred 
package in IRAF, applying a multiplicative gain correc- 
tion to equalize the difference in gain between the two 
readout electronics chains. We then divided by a nor- 
malized flat field exposure to correct for fringing in the 
red-most regions of the spectrum, which has a signifi- 
cant effect on the orders we used for the radial velocity 
analysis. 

At the time of our observations, the TRES detector 
suffered from a very high energetic particle hit rate, the 
source of which is under investigation. The particle hits 
were removed by median-combining each set of multiple 
exposures for the flat fields, target and calibrations, and 
then using the statistics of the median frame to apply 
an upper envelope clip to the target frames themselves. 
This allowed us to avoid having to combine our target ex- 
posures to eliminate the particle hits, which would smear 
out the radial velocity variations slightly as a result of 
the extremely short orbital period of GJ 3236. 

The stability of the TRES instrument allowed us to 
define the aperture trace for each fiber using the high 
signal-to-noise flat fields, which were then used to ex- 
tract the target and calibration spectra, employing the 
"o ptimal" (in a least squares sense) weighting scheme 
of iHewett erall (|1985D . The spectra were wavelength- 
calibrated using the standard ecidentify procedure, 
which employs the (known) dispersion relation of the 
echelle to solve for a wavelength solution over all orders 
of the spectrum simultaneously. This is necessary due to 
the relative paucity of sufficiently intense lines from the 
ThAr lamp in each order of the spectrum at wavelengths 
> 7500 A. 

2.5. Summary of system properties 

Table [S] summarizes the known system properties, from 
the literature (principally the proper motion survey of 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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TABLE 2 

y-BAND LIGHT CURVE OF GJ 3236. 



HJD 


V 


Error'' 




FWHM (pix)': 


Airmass 


X (pix)<l 


y (pix)'^ 


2454734.882776 


14.2717 


0.0030 


-0.137 


6.29 


1.30113 


953.69 


1372.35 


2454734.884327 


14.2676 


0.0029 


-0.094 


6.08 


1.29916 


953.59 


1372.41 


2454734.885866 


14.2672 


0.0029 


-0.068 


5.92 


1.29726 


953.64 


1372.36 


2454734.887406 


14.2698 


0.0029 


-0.112 


6.14 


1.29541 


953.69 


1372.47 


2454734.888957 


14.2696 


0.0029 


-0.093 


6.05 


1.29360 


953.66 


1372.63 



Note. — Tablc[2]is published in its entirety in the electronic edition of the Astrophysical Journal. 
A portion is shown here for guidance regarding its form and content. 

^ Estimated using a standard CCD noise model, including contributions from Poisson noise in the 
stellar counts, sky noise, readout noise and errors in the sky background estimation. Correction 
to the frame magnitude zero-point applied in the differential photometry procedure. More negative 
numbers indicate greater losses. Median FWHM of the stellar images on the frame. The plate scale 
was O'.'67/pix.'^ x and y pixel coordinates of GJ 3236 on the CCD image, derived using a standard 
intensity-weighted moments analysis. 



TABLE 3 

J-BAND LIGHT CURVE OF GJ 3236. 



HJD 


J 


Error^ 


2454513.6413 


10.015 


0.020 


2454513.6421 


10.040 


0.020 


2454513.6430 


10.052 


0.020 


2454513.6438 


10.078 


0.020 


2454513.6447 


10.042 


0.020 



Note. — Table [3] is published in its entirety in 
the electronic edition of the Astrophysical Journal. A 
portion is shown here for guidance regarding its form 
and content. 

^ Equal observational errors of 0.02 mag were as- 
sumed on each data point. 



iLepine fc Sharal 120051 used to select the target stars for 
the MEarth survey) and our own V and MEarth photom- 
etry, converted to the standard Johnson-Cousins system. 

The average color of ~ Ic — 2.73 ± 0.07 indicates 
an average spectral type of approximately M4 using the 
color to spectral type relation of Lcggctt (1992 ) for young 
Galactic disk stars. 

3. INITIAL LIGHT CURVE ANALYSIS (EBOP) 

For detached eclipsing binaries with circular orbits, the 
radial velocity (RV) and light curve models are largely in- 
dependent. We therefore carried out a preliminary anal- 
ysis of the available light curves before starting to obtain 
spectroscopic observations. The principal purpose of do- 
ing this was to determine an extremely precise orbital 
period to better-target the radial velocity observations, 
but we can also constrain the system eccentricity by us- 
ing the phase of the secondary eclipses relative to the 
primary eclipses (related to ecosw). A simple geomet- 
ric argument can then be used to obtain a limit on the 
eccentricity itself, assu ming no a priori knowledge of a;. 

We used jktebop (Southw orth. Maxted k SmallevI 
I2004at ISouthworth et al...2004h) . a mo dified version of 
EBOP (Eclipsing Binary Orbit Program; iPopper &: Etzell 
119811: [Etze3ll980[ ). to obtain these parameters by fitting 
the MEarth light curve. The program was modified to 
fit simultaneously for the EB model, and a simple form 
for the synchronized out-of-cclipse modulations, assum- 
ing they can be approximated by a truncated Fourier 



series: 

m{t) — ai sm{ilt) + bi cos(i7t) + 02 sin(2rit) -|- 62 cos(2il/;) 

(2) 

where fl = ^ir/P, and P is the orbital period of the bi- 
nary. The normalization term is omitted since this is al- 
ready taken into account by the standard EBOP model. 
The revised code yields a fit with very small residuals 
for the present system, and the use of the form in Eq. 

keeps the number of parameters required to describe 
the modulations to a minimum compared to a full spot 
model. Doing so will not necessarily reproduce the cor- 
rect eclipse depths for a given set of physical parameters 
(or vice versa), but should be adequate for deriving an 
accurate system ephemeris. We relax this assumption in 
f|5l where we adopt a full spot model for determining the 
geometric and radiative parameters of the system. 

The following parameters were allowed to vary in the 
fit: the sum of the radii divided by the semimajor axis, 
'"i + 1^2 (where rj — Rj/a, Rj is the stellar radius of 
star J and a is the semimajor axis), orbital inclination 
i, ecosw, central surface brightness ratio J (essentially 
equal to the ratio of the primary and secondary eclipse 
depths), orbital period P, orbital phase zero-point to-, 
light curve normalization, and the parameters ai, 02, 61 
and &2 describing the out of eclipse modulations. Single- 
band light curves showing grazing eclipses constrain the 
ratio of the radii only very weakly (or equivalently, the 
luminosity ratio) , so these parameters and the mass ratio 
(used for computing the oblateness of the stars and the 
reflection effect) were fixed at 1.0 for the initial analysis. 
The results for the orbital parameters were found not to 
change significantly upon revising these to the measured 
values from the spectroscopy once they were available. 

Our light curves are not of sufficient quality to fit 
for the stellar limb darkening, so we assumed a square 
root limb-darkening law with coefficients appropriate to 
the / (Cousins) passband from Claret (2000) using the 
PHOENIX model atmospheres. Given the lack of a spec- 
troscopic constraint on the effective temperature, surface 
gravity and metallicity for either star, we have instead 
assumed solar metallicity, and derived temperatures and 
surface gravities of Ti = 3280 K and log 171 = 4.9, and 
for the secondary, T2 = 3200 K and log 52 = 5.0 by iter- 
ating the modeling process. The assumed temperatures 
were derived using the masses from the combined orbital 
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TABLE 4 

Barycentric radial velocity measurements of GJ 3236. 



HJD 



Absorption (Barnard's star)'' 
vi (km s^^) V2 (km s^-*-) 



Emission (GJ 856A) 
vi (km s^-"-) '02 (km s^-"-) 



Emission (GJ 3379) 
vi (km s^-^) V2 (km s^-"-) 



24547,55 9072 






—44.23 


86.66 


—42.91 


87.65 


2454755.9154 






—49.99 


92.02 


—48.28 


91.13 


2454755.9241 






—53.77 


100.02 


—52.82 


99.68 


2454757 81 1 5 


55.29 


—37.60 


55.70 


—44.31 


56.17 


—42.63 


2454757.8192 


58.96 


—49.14 


61.49 


—51.70 


62.63 


—51.13 


2454757 8272 


61.92 


—56.41 


65.03 


—57.34 


66.39 


—57.10 


2454757.8352 


66.54 


—66.09 


70.44 


—63.70 


71.87 


—62.77 


2454756.8324 






—66.67 


118.68 


—66.39 


119.91 


2454756 8406 






—64.81 


114.28 


—63.84 


116.20 


2454756.8483 






—63.87 


114.44 


—62.96 


117.50 


2454729.8825 






—50.95 


96.13 


—49.89 


96.98 


2454729.8883 


—53.26 


81.21 


—48.21 


90.02 


—47.42 


91.53 


2454729.8941 


—50.49 


79.39 


—44.90 


87.55 


—44.18 


88.84 


2454729 8999 


—46.09 


84.13 


—40.92 


81.58 


—40.06 


83.26 


2454729.9057 


—38.67 


74.80 


—40.44 


77.21 


—39.52 


79.13 


2454759 91 37 






—59.96 


119.56 


—58.57 


123.32 


2454759.9252 






—61.53 


104.53 


—55.20 


102.53 


2454759.9513 






-57.89 


106.49 


-56.91 


107.42 


2454758.8510 


51.55 


-31.67 


48.34 


-37.76 


48.98 


-35.47 


2454758.8587 


48.84 


-26.92 


44.02 


-32.35 


44.54 


-30.38 


2454758.8668 


48.57 


-21.32 


39.19 


-21.43 


40.09 


-18.55 


2454758.8745 


38.38 


-17.90 


33.12 


-16.09 


34.85 


-12.88 


2454730.9731 


90.88 


-101.58 


100.49 


-101.59 


102.23 


-100.68 


2454730.9768 


79.42 


-110.01 


101.24 


-96.71 


102.05 


-96.08 


2454730.9806 


91.46 


-105.07 


99.63 


-93.68 


100.75 


-93.28 


2454730.9843 


90.40 


-110.53 


99.04 


-97.74 


100.43 


-97.25 


2454730.9880 


87.48 


-96.16 


97.66 


-96.19 


99.13 


-95.82 


2454731.9797 


-42.46 


73.24 


-34.37 


74.11 


-33.33 


74.37 


2454731.9835 


-39.77 


73.55 


-36.22 


76.30 


-35.22 


77.22 


2454731.9872 


-46.24 


83.24 


-37.64 


78.70 


-36.75 


80.16 


2454731.9909 


-30.10 


76.54 


-40.50 


83.11 


-39.25 


84.72 


2454732.0024 






-45.76 


91.06 


-44.87 


92.26 



Radial velocities are reported only for epochs where there was a sufficient signal-to- 
noise ratio to obtain a usable cross-correlation. This was satisfied for all the epochs in the 
emission line analysis but only for 21 epochs in the absorption line analysis, due to the 
higher signal-to-noise of the strong Ha lines compared to the continuum. 



and light curve solution, interpolati ng between the com- 
pilation of values for field stars bv lLeggefr3 ()1992D . See 
also S|6] where these temperatures are verified using the 
final model. A gravity darkening exponent oi P — 0.32, 
a valiie app ropriate for stars with convective envelopes 
was also assumed, and the option in jk- 
TEBOP to calculate the reflection effect was used rather 
than fitting for it. 

We report only the orbital parameters (P, to and 
ecosoj) from the ebop analysis, and have adopted these 
for the remainder of this work. 

The parameter uncertainties were derived us- 
ing a Monte Ca.rlo al gorithm built-in to jktebop 
(jSouthworth et al] I2OO50 . Briefly, this algorithm uses 
the best fit to generate a synthetic light curve, injecting 
Gaussian noise with amplitude determined by the 
observational errors (which were scaled such that the 
reduced of the fit was equal to unity) to produce a 
simulated light curve, which is then fit to determine a 
new set of parameters. The parameter uncertainties can 
be estimated using the distribution of the parameters 
from a large num ber of realizations (here, 10 000) of this 
process. See also ISouthworth et all (|2004allbl ) for more 
details. 

Times of minimum light derived for the 11 well- 
sampled eclipses in our V and MEarth observations are 
reported in Table [71 These were estimated using the 



method of lKwee fc van WoerdenI (|1956( l over a ±0.03 re- 
gion in normalized orbital phase around each eclipse. We 
subtracted the truncated Fourier series in Eq. ([2]) from 
the light curves before applying this analysis to reduce 
the effect of any systematic errors arising from the out- 
of-eclipse modulation. 

4. RADIAL VELOCITY ANALYSIS 

Radial velocities were obtained using the two- 
dimensional cross-cor relation algorithm todcor 
(jZucker fc Mazel]|1994f ). which uses templates matched 
to each component of a spectroscopic binary to si- 
multaneously derive the velocities of both stars, and 
importantly for the present application, an estimate of 
their luminosity ratio (a) in the spectral bandpass. 

We performed both a standard absorption line based 
cross-correlation analysis, and an emission line analysis 
using the strong Ha emission observed in both compo- 
nents of GJ 3236. The signal to noise ratio of the lat- 
ter is substantially higher in our data, but we caution 
that the source of the Ha emission in our target stars is 
not well-understood. The emission line analysis results 
could suffer systematic errors depending on the relative 
velocities of the regions of the chromosphere emitting 
the Ha and the stellar photospheres, and this is exac- 
erbated in the present system because it is highly likely 
that the stellar spin and binary orbit are synchronized, 
both from the light curves and expectations from tidal 
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Fig. 1. — Images of GJ 3236 centered on its position as measured 
from the MEarth data. The circle shows the approximate position 
and size of the 5" (radius) photometric aperture used to derive our 
hght curves. Data are from the first and second epoch Palomar sky 
surveys as provided by the Digitized Sky Survey (top and center 
panels), and the MEarth stacked master image (bottom panel). 
The approximate epochs of the images are 1954.1 (POSS-1), 1992.7 
(POSS-2), and 2008.8 (MEarth). All three panels have the same 
center, scale and alignment on-sky, with north up and east to the 
left, covering 1' in the horizontal direction. 

theory (jZahnl [19771 ) . We therefore conservatively adopt 
the results from the absorption line analysis despite the 
larger errors. 

The absorption line analysis used as a template 
spectrum a single epoch observation of Barnard's star 
(Barnard 1916; also known as GJ 699) taken on UT 2008 
October 20. We used a wavelength range of 8700 to 8850 
A in order 50 of the spectrum to derive the velocities, 
since this region contains a number of reasonably strong 
metallic lines in M-dwarfs, and is free of telluric absorp- 
tion lines. The Ca II infra-red triplet at 8498, 8542, and 
8662 A was deliberately avoided since the absorption in 
these lines is saturated, and they are therefore intrinsi- 
cally very broad, and emission cores are often observed in 



TABLE 5 

Summary of the photometric and astrometric 
properties of the g,j 3236 system. 



Parameter 


Value 


a 1-1 
"2000 ' 


03''37"'14''.08 




-|-69°10'49'.'8 


fj,a COS 


0'.'142 yr-l 




-0'.'132 yr-l 


MEarth'^ 


11.05 ±0.05 




14.28 ±0.05 




11.55 ±0.05 




9.806 ± 0.023 


-H'2MASS'^ 


9.198 ± 0.028 


.?f2MASS'* 


8.967 ±0.022 



Equinox J 2000 0, epoch 2000.0.*= From 
iLepine fc Sharal ()2005l ).'^ Median of the measured 
out-of-eclipse fluxes from our light curves. We esti- 
mate that the observational errors are dominated by 
the uncertainty in the photometric calibration, which 
is approximately 0.05 mag for both passbands. Ic was 
computed from the measured MEarth and Vj magni- 
tudes using Eq. ^■'^ We quote the combined un- 
certainties from the 2MASS catalog, noting that the 
intrinsic variability of our target means in practice 
that these are underestimates. 

these lines in active stars, which would severely compli- 
cate the use of an inactive field star as a template. We as- 
sumed a barycentric radial velocity of — 110.13 km s^^ for 
Barnard's star, derived from presently unpublished CfA 
Digital Speedometer measurements spanning 17 years. 

For the emission line analysis, we used as templates 
spectra of two M4 dwarfs that were found to display Ha 
emission, GJ 856A and GJ 3379, observed as part of 
another program using the same instrument. GJ 856 A 
was observed on UT 2008 October 17 and GJ 3379 was 
observed on UT 2008 October 19. We used a wavelength 
range of 6520 to 6645 A in order 37 for the emission line 
cross-correlations, noting that the cross-correlations are 
dominated by the Ha lines. 

Barycentric radial velocities of GJ 856A and GJ 3379 
were determined by cross-correlation with Barnard's 
star, using the 8700 to 8850 A region as for our absorp- 
tion line analysis of GJ 3236. 

Given that todcor provides a measure of the corre- 
lation between the observed and template spectra, it is 
possible to run cross-correlations against a series of tem- 
plates and use the correlation value as an indicator of the 
template match. Both the effective temperature and the 
rotational velocity of the templates affect the solution, 
so it is vital to explore this parameter space. Due to 
the availability of a limited range of observed templates, 
the effective temperature could not be varied, but we 
were able to rotationally broaden the template spectra 
in order to create primary and secondary templates with 
Wrot G {1,2,4,6,8,10,12,16,20,25,30,35,40} km s^^. 
We used linear limb darkeni ng coe fficients ui — 0.5782 
and U2 = 0.5934 from Clarell pOOl for the SDSS z pass- 
band (a close match to the effective wavelength of the 
spectroscopy), which correspond to our assumed tem- 
peratures for the primary and secondary (see 
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TABLE 6 

Derived orbital parameters for, the GJ 3236 system. 



TABLE 7 

Measured times of minimum light for GJ 3236. 







Value^ 


P 

to 




0.7712600 ± 0.0000023 days 
2454734.99586 ± 0.00012 HJD*' 


e cos uj 
e 




0.00008 ± 0.00020 

< 0.0022 (95% CL)<: 

< 0.0112 (99% CL)'= 


a 




0.60 ± 0.04 


Absorption 


line solution (adopted) 


7 




10.06 ±0.94 km s"! 
85.6 ±2.1 km s"! 
114.8 ± 1.9 km s-l 


q 

asini 
Ml sin^ i 
M2 sin^ i 




0.746 ± 0.023 
3.053 ± 0.044 Rq 
0.368 ±0.015 M0 
0.275 ± 0.014 M0 






6.3 km s-i<i 
5.6 km s"-*^ 


Ha emission line solution"^: 


GJ 856A 


7 

K2 




12.87 ±0.19 km s^^ 
88.48 ±0.33 km s"! 
114.71 ±0.49 km s"! 


1 

a sin i 
Ml sin^ i 
M2 sin^ i 




0.7713 ± 0.0045 
3.0962 ± 0.0091 R0 
0.3786 ± 0.0037 Mq 
0.2919 ±0.0025 Mq 






1.3 km s-i 
1.9 km s-l 


Hq emission line solution"^: 


GJ 3379 


7 

K2 




13.98 ±0.20 km s^^ 
88.65 ±0.34 km s"! 
114.87 ±0.52 km s"! 


1 

a sin i 
Ml sin^ i 
M2 sin^ i 




0.7717 ± 0.0047 
3.1011 ±0.0096 R0 
0.3802 ± 0.0039 Mq 
0.2934 ± 0.0026 Mq 


o"l 
0-2 




1.3 km s~^ 
2.0 km s-l 



^ We report 68.3% confidence intervals, with error 
bars symmetrized by adopting the mean of the 15.85% 
and 85.15% levels as the central value.'' Ephemeris 
zero point, chosen to correspond to the epoch of 
the first primary eclipse in the T^-band data from 
^12. 21 " Derived assuming only the measured ecosw 
and a uniform distribution in w, for 95% and 99% 
confidence levels. These confidence levels are in fact 
lower limits, since there were no noticeable differences 
in the eclipse durations, or deviations from a circular 
orbit in radial velocity, e = was assumed hence- 
forth for the radial velocity modeling.'* RMS of the 
residuals after subtracting the Keplerian orbit model 
from the data. These are representative of the typical 
uncertainty on each RV point." We re-emphasize that 
the Ha radial velocities, despite having smaller ran- 
dom errors than the absorption line velocities, may 
have an unknown and potentially significant system- 
atic error, due to the uncertainty in the source of 
the Ha emission line, and whether it tracks the stel- 
lar photosphere. We therefore conservatively adopt 
the absorption line solution for the remainder of the 



HJD 


(O - C) (s) 


Cycle^ 


Band 


2454734.995569 


-32.5 ±52.0 


0.0 


V 


2454736.923700 


-34.5 ±26.7 


2.5 


V 


2454741.937128 


-7.8 ±17.0 


9.0 


MEarth 


2454742.708283 


-19.6 ± 18.9 


10.0 


MEarth 


2454745.793562 


7.2 ±18.9 


14.0 


MEarth 


2454762.761020 


-22.0 ±25.2 


36.0 


MEarth 


2454765.846327 


7.9 ±78.0 


40.0 


MEarth 


2454766.617531 


1.7 ± 13.7 


41.0 


MEarth 


2454767.774152 


-28.4 ±67.0 


42.5 


MEarth 


2454770.859614 


18.9 ±36.6 


46.5 


MEarth 


2454771.630638 


-7.6 ±38.6 


47.5 


MEarth 



^ Cycle number, counting from at the primary 
eclipse at t^, in units of the orbital period. Integer 
numbers correspond to primary eclipses. 



We ran todcor on the grid of 169 template combina- 
tions, and sought the maximum of the resulting correla- 
tion values to determine the most likely Wrot sin i values 
for the two components of the binary system. Barnard's 
star is a very slow rotator, so the rotation of the template 
is negligible given our spectral resolution. Our analy- 
sis indicates frot,i sinz = 25 km s^^, which is consistent 
with the expected value of Wrot,i sinz = 24.6 km s^^ if the 
spin and orbit are synchronized. Wrot,2sini was poorly- 
constrained, with the maximum correlation occurring at 
3 km s~^, whereas we would expect 19.5 km s~^ assum- 
ing synchronization. We have simply adopted the latter 
value, since the light curves and primary velocity indi- 
cate that the system is indeed synchronized. Using the 
templates with iirot.isini — 25 km s^^ and Wi.ot,2sini — 
20 km s~\ TODCOR indicates a = 0.60 ± 0.04. 

For the emission line analysis, we assumed the stellar 
rotation values from the absorption line analysis, rota- 
tionally broadening the observed spectra of GJ 856A and 
GJ 3379. The intrinsic rotational broadening in both of 
these templates is again negligible for our purposes. TOD- 
COR gave a maximum correlation for ana — 0.52 ± 0.02, 
which provides a measure of the relative Ha emission 
luminosities of the two stars. 

Since the orbital period and time of primary eclipse are 
extremely well-determined from the light curves, we held 
these values fixed in the final radial velocity orbital solu- 
tion. We also assumed a circular orbit. Figure [2] shows 
the resulting radial velocity curves. The parameters de- 
rived from the radial velocity analysis are given in Table 
[6l We assumed equal observational errors on each radial 
velocity point, scaled such that the reduced of the fit 
was equal to unity. 

Radial velocities of the components of double-lined 
spectroscopic binaries are prone to a "peak pulling" ef- 
fect as they approach the 7 velocity (at times of mini- 
mum separation of the two components in the spectrum) . 
The TODCOR method can largely eliminate this, provided 
the template spectra match the target sufhciently well. 
In our experience, this effect is particularly sensitive to 
errors in the assumed rotational broadening: if the tem- 
plate contains too much or too little rotation, the veloci- 
ties would be pushed toward or away from 7 to compen- 
sate. Examining the residuals of our solution in Figure O 
indicates little or no "peak pulling" , which further justi- 
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Fig. 2. — Phase-folded radial velocity curves for the GJ 3236 sys- 
tem. Shown in each panel are the radial velocities for the primary 
(filled symbols) and the secondary (open symbols), with the best- 
fitting Keplerian orbit overplotted using solid and dashed lines for 
the two components, respectively. The lower sub-panels show the 
residuals after subtracting the model from the data. Top: absorp- 
tion line solution cross-correlating against Barnard's star; Center: 
emission line solution using GJ 856A; Bottom: emission line solu- 
tion using GJ 3379 as the template. 

fies our assumed rotational velocities and the assumption 
of negligible rotation in the templates. 

5. LIGHT CURVE ANALYSIS (WILSON & DEVINNEY 
CODE) 

In order to properly account for spots on the com- 
ponents of GJ 3236, we u se the most recent v ersion 
(WD2007) of the standard IWilson fc DevinnevI (|197lD 
code to derive the geometric and radiative parameters of 



the system. This code incorporates a full physical model 
including spots and usin g model atmospheres. 

WD normally uses the lKurucj (|1979l . |1993[ ) model at- 
mospheres to calculate the emergent flux from surface 
elements on each binary component. These models span 
an effective temperature range of 3500 — 50000 K, and for 
temperatures below 3500 K, WD interpolates between 
a black body of the appropriate temperature and the 
3500 K model atmosphere. This is weighted in such a 
fashion as to produce a pure black body at 1500 K, and 
for intermediate temperatures there is a smooth ramp 
toward a pure model atmosphere at 3500 K. Molecu- 
lar sources of opacity increasingly begin to dominate the 
optical spectra of M-dwarfs in this temperature range, 
so a black body is clearly a poor approximation for the 
emergent flux. Since both components of GJ 3236 ap- 
pear to be cooler than 3500 K, we have modifled our 
copy of W D to substitute phoenix NextGen model at- 
mospheres ([Hauschildt. AUard &: BaronI I1999D in place 
of these interpolations for effective temperatures cooler 
than 3500 K. 

5.1. Assumptions 

Our assumptions for the WD light curve model are 
detailed in Table [51 We have assumed zero orbital eccen- 
tricity, synchronized spin and orbit for both components 
of the binary, and no third light, which are justified in 
gOand sp 

The limited phase coverage and lack of repeat obser- 
vations in our y-band photometry (leading to an un- 
certainty in any corrections for color-dependent "second- 
order" atmospheric extinction), and poor precision (large 
scatter) in the J-band photometry means in practice that 
these are not useful in aiding the determination of the 
geometric (and thus physical) parameters of the system. 
Our experiments with the WD model indicate that these 
provide minimal additional constraints on any of the 
parameters when compared to using the MEarth data 
alone. We have therefore elected to use only the MEarth 
data for fitting for the system parameters, and then eval- 
uated this model with respect to the other passbands as 
a check. 

Doing so necessitates fixing one of the component effec- 
tive temperatures, since the relative eclipse depths mea- 
sured in a single passband merely determine the relative 
temperatures of the two components. We have there- 
fore assumed Ti — 3280 K following SjS] In practice even 
given multi-band photometry it is usually necessary to 
fix this parameter: although the color-dependence of the 
eclipse depths does in principle provide information on 
the value of Ti, this is normally poorly constrained due 
to the difficulty of measuring differential colors to high 
precision. 

In order to reproduce the out-of-eclipse variations seen 
in our optical light curves, at least two spots must be 
introduced into the model. The configuration of these 
spots affects the measured eclipse depths, depending on 
whether a spotted or unspotted part of the photosphere 
is eclipsed. This predominantly affects the derived ratio 
of surface brightnesses as measured by the ratio of depths 
of the primary and secondary eclipses (and thus the ra- 
tios of the component radii and effective temperatures 
inferred from this), and a combination of (i?i -I- i?2)/a 
and the orbital inclination inferred from the total eclipse 
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TABLE 8 

Parameters used in the WD light curve model. 



Parameter 


WD name 


Value or "varied" 


Description 


Binary type 




detached 




to 


HJDO 


2454734.99586 (HJD) 


Phase zero-point^ 


P 


PERIOD 


0.77126 days 


Orbital period^ 


dP/dt 


DPDT 





First derivative of period 


a, sin i 


SMA X sin i 


3.053 ± 0.044 Rq 


Projected semimajor axis^'^ 


n 
H 


RM 


0.746 ± 0.023 


Mass ratio^ 


T 


VGA 


10.06 ±0.94 km s"! 


Systemic velocity'' 






Q 






Fl 


1.0 


Primary rotation parameter 




F2 


1.0 


Secondary rotation parameter 




HLA 


varied 


Light curve normalization 


Tx 


TAVH 


3280 K 


Primary effective temperature 


To 


TAVC 


V CJbL l"\A 




Qi 


PHSV 


varied 


Primary surface potential"^ 


Ha 


PCSV 


varied 


Secondary surface potential'^ 


i 


INCL 


varied 


Orbital inclination 




ELS 





Third lio-ht 


Ai 


ALBl 


0.5 


Primary bolometric albedo'^ 


Ai 


ALB2 


0.5 


Secondary bolometric albedo'* 


I3i 


GRl 


0.32 


Primary gravity darkening exponenf* 


ft 


GR2 


0.32 


Secondary gravity darkening exponenf* 


0s,l 


XL ATI 


60° 


Spot 1 latitude" 


4>s,i 


XLONGl 


varied 


Spot 1 longitude 


rs,i 


RADSPl 


varied 


Spot 1 radius 


Ts,i/Ti 


TEMSPl 


0.96 


Spot 1 temperature contrast 


es,2 


XLAT2 


60° 


Spot 2 latitude'' 




XL0NG2 


varied 


Spot 2 longitude 


^3,2 


RADSP2 


varied 


Spot 2 radius 




TEMSP2 


0.96 


Spot 2 temperature contrast 


a 




0.60 ±0.04 


Light ratio'' 



From Table Although the WD parameter SMA is the semimajor axis itself, 
we fixed asini in our MCMC procedure to the value from the radial velocities, 
calculating t he required SMA value using i (iNCL) at each iteration. Modified 
iKopall (119591 ) potentials. These are roughly inversely proportional to stellar radii 
for this detached binary; for a clear description of their definition, we refer the 
reader to the documen tation for the new WD graphical user interface phoebe 
()Prsa fc Zwitteii [20051 ).'* Values appropriate for convectivc atmospheres. See 
text. 



depth. 

Since our V^-band data do not cover the portions of 
the light curve out of eclipse, and the scatter in the J- 
band data is too large to usefully constrain the model, 
the spot modeling must be done using single-band light 
curve information. Spots are introduced in the WD code 
using a standard four-parameter model, where each spot 
is parametrized using the spot latitude 9s and longitude 
4>s, radius and temperature contrast Tg/Tp. The latter 
quantity gives the ratio of effective temperatures of the 
spotted and unspotted photosphere. For a single-band 
light curve, the spot longitude is well-constrained by the 
phase of the out of eclipse modulation, but the combi- 
nation of 0s, Tg and Tg/Tp is largely degenerate (these 
quantities determine the amplitude and shape of the vari- 
ation). Moreover, we can place each spot on either star. 
Although the data do rule out some of the possible spot 
locations, there are still a range of possibilities, which 



introduces a systematic error into the geometric param- 
eters of the system. In principle, measuring the color- 
dependence of the out-of-eclipse modulation would allow 
this degeneracy to be reduced. 

Given the presently available data, we instead probe 
this systematic error by considering three characteris- 
tic spot configurations that are consistent with the light 
curve. These are: (1) one close to polar spot on each star, 
with 9s — 60°, with both spots located in the same hemi- 
sphere as the chord traversed during the eclipse, result- 
ing in the spot on the primary being eclipsed by the sec- 
ondary during primary eclipse, and vice versa during sec- 
ondary eclipse; (2) likewise, only with the spots located 
in the opposite hemisphere, resulting in no spots being 
eclipsed; and (3) a case intermediate between these, with 
a "symmetric" spot configuration of two spots on each 
star, at latitudes of 60° with one located in each hemi- 
sphere. In each model, the spot latitudes and tempera- 
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ture ratios were held fixed, allowing only longitude and 
radius to vary, and in configuration (3) we enforced sym- 
metry of the spot pattern about the equator, i.e. both 
spots on one star were forced to have the same longi- 
tudes and radii. This means that all three models have 
the same number of parameters (4). 

5.2. Model fitting and error estimates 

To derive robust error estimates, including the cor- 
relations between parameters, we have implemented 
a Markov Chain Monte Carlo (MCMC) algo rithm 
aroun d the WD light curve model. Following iFordI 
2005f). we used the Metropolis-Hastings algorithm 



Metropolis et al.| [l953: Hastings 1970) to estimate the 



a posteriori joint probability distribution of the fitted 
model parameters. The fitting statistic was the tradi- 
tional J where the observational error estimates were 
scaled such that each light curve's reduced was equal 
to unity. 

We briefly summarize the Metropolis-Hastings algo- 
rithm here. Starting from an initial point in parameter 
space, the algorithm takes the most recent set of param- 
eters and perturbs one or more parameters by a random 
Gaussian deviate. If the perturbed parameter set has 
a lower than its progenitor, it is accepted as a new 
point in the chain. If it has a larger it is accepted 
with a probability exp{— Ax^/2). If it is not accepted, 
the original point is repeated in the chain. The size of 
the perturbations were adjusted so that 20 — 30% of the 
proposed points were accepted. 

We included the spectroscopic light ratio between the 
primary and secondary components as an observation in 
our x^ statistic, which in practice was implemented by 
adding an extra contribution to of the form: 



Apnor 



^measured — CtWD 



(3) 



where ameasured and (Tq, are the observed light ratio and 
its error, taken from Table [51 and awD is the value calcu- 
lated from the WD model. This is necessary for systems 
with near-circular orbits and grazing eclipses because the 
light curves only very weakly constrain this quantity (or 
equivalently, the ratio of the component radii). Uninfor- 
mative (uniform) priors were assumed on all other pa- 
rameters. 

The final parameters and error estimates were deter- 
mined by taking the 68.3% confidence interval (15.85 and 
84.15 percentiles, corresponding to ±1 standard devia- 
tion for a Gaussian distribution) of the samples gener- 
ated by our MCMC procedure. The chains were run un- 
til they contained at least 10^ samples, and we discarded 
the first 10% of the points in each chain in order to en- 
sure they had converged. These parameters are reported 
in Table [HI and Figure [3] shows the light curves with the 
best-fitting model overplotted. 

5.3. Discussion 

The most serious limitations in the present analysis 
result from the use of only single-band light curves to 
model the out of eclipse variations. This is clear from 
comparing the results for the three solutions reported 
in Table [51 where wc find significant (> 1 standard de- 
viation) differences in the orbital inclination i, ratio of 



1,20 
1-25 
1.3D 
1.35 
-0.03 
-0.02 
-0.01 

o.oc 



14.25 
U.JO 
14.35 
14.40 
14.45 
14.50 
-0.0^ 
-0.03 
-0.02 
-0.01 
0.00 
0.01 
0.02 









- 


- 

y 


























^ 




'A^r^^ i^^A^ 










V 















Fig. 3. — Phase-folded light curves in MEarth, V and J pass- 
bands (top, middle and bottom panels, respectively). In each 
panel, the upper sub-panels show the light curve (black points) 
with the WD model 3 (see Table [9)l overplotted (solid lines), and 
the lower sub-panels show the residuals after subtracting the model 
from the data. The uppermost panel for the MEarth data shows a 
binned light curve (in 250 bins each spanning 0.004 in normalized 
orbital phase) to allow the features to be more clearly distinguished 
given the large number of data points. 



effective temperatures T2/T1, and in (i?i -f i?2)/a, which 
depend on the assumptions made regarding the locations 
of the spots on the components of the binary. 

Given our radial velocity errors, the orbital inclina- 
tion uncertainty has little effect on the final physical pa- 
rameters, as evidenced by the component masses and 
semimajor axis reported in Table[9l which are essentially 
identical to within the observational errors. The domi- 
nant uncertainty in the final solution is therefore in the 
component radii. 

Of the three models considered, the symmetric spot 
distribution in model (3) seems the most physically 
realistic given the known properties of spots on low- 
mass stars. These are often found to be polar in 
light curve models (e.g. iRodono et al.' '1986'). More- 
over, Doppler imaging studies ( e.g~[Barncs et al. 200(|; 
[Barnes fc Collier Cameronll2001f ) indicate that in reality, 
the surfaces of low-mass stars probably have many small 
spots distributed over a range of latitudes. This situa- 
tion is largely indistinguishable from a small number of 
large, polar spots in one-dimensional light curve data, 
although it is important to note that the effects on the 
measured eclipse depths could be different. There is no 
a-priori reason to expect the spots to be concentrated 
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TABLE 9 

Derived geometric, radiative and physical parameters of the GJ 3236 system. 



Parameter 


Model 1 


Value^ 
Model 2 


Model 3 


HLA 




U.UUUoUo/^ ± U.UUUUU14 


c\ nnnom c; _l n c\r\c\c\(\c^ci. 

u.uuuzyio ± u.uuuuudd 


T2/T1 


0.9801 ±0.0016 


0.97914 ± 0.00085 


0.9774 ±0.0014 


^1 


9.01 ±0.12 


8.781 ±0.030 


8.794 ±0.098 




0.39 ± U.14 


9.23d ± 0.041 


8.806 ± 0.156 


i 


OZ.OUO ± U.U4i 


QQ "701 _1_ n OQO 

83. (21 ± 0.032 


83.154 ± 0.U4/ 


<t>s,l 


-17.2° ± 1.5° 


-14.0° ± 1.4° 


-21.6° ± 2.0° 


rs,l 


36.5° ± 1.1° 


43.83° ±0.41° 


25.90° ± 0.42° 


<Ps,2 


1 T c: _L 1 AO 


on AO _L 1 00 


00 no 1 1 10 
ZZ.U lb i.l 


rs,2 


4d.7 ± 1.0 


59.95 ± U.9o 


34.54 ± U.d5 


(Hi + H2j/a 


n ooQiT? _L n nnnoo 
U.iZoy/ ± U.UUUoz 


n 01 airx _L n nnn'zQ 
U.ZiOfU ± U.UUU(0 


n 001 Q^z _L n nnnoc; 
U.ZZLOi ± U.UUUo5 


R2/R1 


0.850 ±0.028 


0.7380 ±0.0048 


0.783 ± 0.024 


Ml 


n 'i77 -L n ni 

U.Oi I ^ u.u±u 






M2 


0.282 ± 0.015 


0.280 ± 0.015 


0.281 ± 0.015 Mq 


Ri 


0.3729 ± 0.0078 


0.3829 ± 0.0057 


0.3828 ± 0.0072 Rq 


R2 


0.3167 ±0.0075 


0.2828 ±0.0043 


0.2992 ± 0.0075 R© 


a 


3.077 ± 0.044 


3.071 ± 0.044 


3.075 ± 0.044 R0 


log 51 


4.872 ± 0.015 


4.8456 ± 0.0081 


4.850 ±0.014 


log 52 


4.886 ± 0.019 


4.9819 ±0.0115 


4.935 ± 0.019 


frot,l 


24.44 ± 0.52 


25.12 ±0.37 


25.05 ± 0.50 


frot,2 


20.77 ±0.47 


18.54 ±0.27 


19.60 ±0.44 


-^samp 


1.61 X 10« 


2.35 X 10^ 


1.37 X 106 


""MEarth 


0.0061 


0.0061 


0.0061 




0.0074 


0.0076 


0.0066 




0.027 


0.028 


0.027 



We report 68.3% confidence intervals, with error bars symmetrized 
by adopting the mean of the 15.85% and 85.15% levels as the central 
value.'' RMS scatter of the residuals after subtracting the model from the 
data. 



into one hemisphere given the close to edge-on inclina- 
tion of the binary orbit. We therefore favor model (3), 
but report all three solutions to provide a realistic evalu- 
ation of the systematic error in our results, and reiterate 
that the data do not distinguish between the three pos- 
sibilities, with all of these having comparable values 
and residuals. 

As a check, we compare the model fit to the MEarth 
data alone with our V and J-band data in Figure [3l In 
order to do this, we have refit the model for these bands 
allowing only the light curve normalization parameter 
(HLA) to vary, with all other parameters fixed to the val- 
ues determined from the MEarth data. The model is 
consistent with the J-band data within the scatter, and 
is reasonably consistent with the l^-band data, consid- 
ering that the early parts of the secondary eclipse curve 
were taken during quite non-photometric conditions, and 
that we have not accounted for color-dependent or differ- 
ential atmospheric extinction effects, which are expected 
to be larger in this band than the redder MEarth or J 
bandpasses, and cannot readily be determined due to the 
lack of repeat observations. 

The presence of spots and proximity effects in very 
close binary systems have a small influence on the shape 
of the radial velocity curve, causing it to depart from 
a Keplerian orbit as we assumed in Sj4l We have in- 
vestigated the influence of these effects on our results 
by comparing radial velocity curves generated using the 
WD model and the simple Keplerian model, finding that 



the corrections to radial velocities taken out of eclipse 
are dominated by the spots, the properties of which 
are largely unknown. However these are very small 
(< 0.1 km s~^) compared to the uncertainties in our 
radial velocity measurements (approximately 6 km s^^), 
so we are therefore justified in neglecting them for the 
present analysis. 

6. EFFECTIVE TEMPERATURES, LUMINOSITIES AND 
KINEMATICS 

Provided the effective temperature of one component 
of an eclipsing binary is known, the combined light curve 
and radial velocity parameters are then sufhcient to de- 
termine intrinsic, bolometric luminosities of both com- 
ponents of the system. This follows directly from the 
definition of the effective temperature in terms of the 
Stefan-Boltzmann law: 

Lboi = 47ri?Vr4 (4) 

By using bolometric corrections and the measured sys- 
tem magnitudes, we can then infer the distance, provided 
the reddening can be constrained or assumed to be zero. 
In the present case, the latter is a re asonable assump- 
tion since our target is very nearby ('G liese &: Jahreisa 
199ll 



nvea "photometric parallax" of 47 ± 8 pc, and 



Lepinel[2005l gives 21.3 ± 4.4 pc; note that both of these 
assume the system is a single star, which means they un- 
derpredict the distance to a near equal luminosity binary 
such as the present one by a factor of approximately v^) . 

In order to determine the effective temperatures, we 
assume the intrinsic colors and bolometric corrections 
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for M-dwarfs compil ed by iBesselll (|1995l . Il991h and 
iBessell fc BrettI ()1988l ) to convert the integrated system 
V — I color presented in Table [5] into Ti , assuming the 
measured effective temperature ratio and radii for both 
components from the light curve mo del. We used the 
polynomial fits of Hillenbrand (jl997f ) , which cover the 
required spectral type range, and assume an uncertainty 
of 0.1 mag in these fits as stated in her Appendix C. We 
also assume a ±100 K systematic uncertainty in the ef- 
fective temperature scale, which has been included and 
propagated in the errors we report. Table [TOl gives our 
derived parameters for the GJ 3236 system. 

The refined value of Ti reported here differs by ap- 
proximately 35 K from the assumed value in used to 
determine limb darkening coefficients. This discrepancy 
is < la, and should have a negligible effect on the pa- 
rameters determined from the light curve analysis, so we 
elected not to repeat this using the updated value of Ti 
given the computational expense involved. 

Given the full system kinematic information (posi- 
tion, proper motions, 7 velocity and parallax from the 
EB analysis) we can also calculate the UVW compo- 
nents of the s pace velocity relative to the S un. We use 
the method of iJohnson fc SoderblomI (|1987f) . but adopt 
the definition that positive U values are away from the 
Galactic center (a left-handed coordinate system) for 
better consistency with the literature. These are in- 
cluded in T able [TUl Compa r ing to the velocity ellipsoids 
derived by iGhiba fc Beeri (|2000( ). GJ 3236 appears to 
be a member of the Galactic disc, and lies within the 
bounds of t he you ng Galactic disc population as defined 
bv iLeggetd (jl99l . Comparing to the UVW space mo- 
tions of various solar neigi iborhood moving groups (e.g. 
ISoderblom fc Mavodll993l ). the velocities are consistent 
with with those of the Hyades moving group, which 
has (U,V,W) = {+2,^ ± 6, -17 ± 6,-11 ± 12) km s^^ 
(jZhao et al.ll2009[ l. Despite this possible membership in 
the Hyades group, we note that the dispersion in stel- 
lar parameters for the group members does not provide 
a very useful con straint on age_or metallicity for the 
present system. Zhao et al. (2009) obtain a metallicity 
of [Fe/H] = -0.09 ± 0.17 for the Hyades group. 

7. DISCUSSION 

Figures m and [5] show the position of GJ 3236 relative 
to other known objects and theoretical models of the 
M-dwarf mass-radius and mass-effective temperature re- 
lationships. 

Although many of the best-observed systems are found 
to exhibit larger radii than the models predict, the slope 
of the mass-radius relation between the components of 
each system (or equivalently, the ratio of the component 
radii) is generally found to agree well between the ob- 
servations and the theoretical models. Examining the 
positions of our solutions for GJ 3236, our models 1 and 
3 for the out of eclipse modulations are reasonably consis- 
tent with this expectation, whereas the slope defined by 
the components for solution 2 is substantially different. 
Although this indicates the assumptions used in model 
2 may be unrealistic, the present observational data for 
GJ 3236 do not distinguish between the three possibili- 
ties. Using this as a constraint in the modeling would be 
dangerous because the results would then no longer yield 
a completely independent test of the theoretical models. 




0,2 0.3 0.4 

M / Mg 

Fig. 4. — Mass-radius relation for eclipsing double-lined spec- 
troscopic binary systems with one or more components below 
0.4 M0. Black points with error bars show the two components 
of the GJ 3236 system, for all three solutions reported in Ta- 
ble [9I The gray points are known systems from the literature 
IIMoTales et al.ll200a IRiba}j |2003 : IVaccaro et"al|[2007l : IBlake et all 
120081) . The black line shows the theoretical mass-radius relation- 
ship from .Baraffc et al., (19981 ) for an age of 1 Gyr. 




Fig. 5. — As Figure |4] but showing the mass-effective tem- 
perature relationship. Again, all three solutions are shown 
for GJ 3236, but the differences are relatively minor in this 
plane. We also include several more m assive M-dwarf EBs in 
this figure for comparison: YY Gem l ITorres &: Ribaj [2002 ), 
GU Boo (LoDcz-Morales & Ribas 20Q1), and OGLE BW3 V38 
(|Maceroni & Montalban ,2004) ■ The black line shows the theo- 
retical relationship from iBaraffe et al.l 1(1993 ) for an age of 1 Gyr. 

so instead we suggest that a better solution to the issue 
would be obtaining improved multi-band photometry to 
better-constrain the spot parameters. 

Therefore, within the present observational uncertain- 
ties, the radii of GJ 3236 appear to be consistent with 
the theoretical predictions at approximately Icr, when we 
account for the systematic uncertainties represented by 
the three solutions we have presented. Furthermore, the 
central values appear to lie above the theoretical curve. 
This is largely consistent with the results for the well- 
known systems in the literature, where the radii are typ- 
ically found to be 10 — 15% larger than the theoretical 
predictions. 

In the effective temperature versus mass plane (Fig- 
ure [5]) , the effective temperatures of the components of 
GJ 3236 are found to be cooler than the models pre- 
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TABLE 10 

Effective temperatures, luminosities and kinematics for GJ 3236. 



Parameter 




Value 






Model 1 


Model 2 


Model 3 


-'I 


oo±o in j-iu 


QQ1 n -t- 1 1 n K 




To 


3247 + 108 K 


3241 ± 108 K 


3238 + 108 K 


-^bol.l 


0.0152 ± 0.0021 L0 


0.0160 ±0.0022 L0 


0.0160 + 0.0021 L0 


-^bol 2 


0.0101 ± 0.0014 L© 


0.0080 ± 0.0011 L0 


0.0089 ± 0.0012 L0 


Mv 


11.17 ±0.30 


11.23 ± 0.30 


11.19 ±0.30 


Mi 


8.44 ±0.27 


8.50 ±0.27 


8.46 ± 0.27 


(m - M) 


3.11 ±0.28 


3.05 ±0.28 


3.09 ± 0.28 


d 


42.2 ± 5.5 pc 


41.1 ±5.3 pc 


41.8 ±5.4 pc 


U 


+34.2 ± 3.9 km s"! 


+33.5 ±3.8 km s"! 


+34.0 ± 3.9 km s~i 


V 


-20.7 ± 4.0 km s"! 


-20.0 ±3.9 km s"! 


-20.5 + 4.0 km s"! 


W 


-2.5 ± 1.7 km 


-2.4 ± 1.7 km 


-2.5 ± 1.7 km s-l 



diet, by approximately 2a. This is in good agreement 
with many of the well-known systems in the literature as 
shown in the figure, such as CM Dra. 

Compared to the well-known systems CM Dra and 
CU Cnc, GJ 3236 has a shorter orbital period, and the 
Ha emission in the spectroscopy and out-of-eclipse mod- 
ulations we observe in the light curve are indicative of 
high activity levels, with both the photometric period 
and the spectroscopic line broadening apparently con- 
sistent with the stellar spin being synchronized t o the 
binar y orbit, as expected from tidal theory (e.g. IZahnI 
fl97l. 

iChabr ler "eTall (pOOTl ) propose two hypotheses to ex- 
plain the observed radius discrepancies between theo- 
retical models and eclipsing binaries: (1) that the in- 
flated radii result from reduced convective efficiency due 
to high rotation rates and large magnetic fields; or (2) 
that magnetic spot coverage of the surfaces leads to re- 
duced heat flux, and thus larger radii and cooler effec- 
tive temperatures. Hypothesis (2) has been found to be 
consistent with some of the well-known e clipsing binary 
systems, e.g. YY Gem (Ivan "Gengil92l iJov fc Sanfo^ 
119261: [Torres fc Ribai |2OO0), where the radius discrep- 
ancy (and the discrepancy in effective temperature) can 
be explained by the presence of starspots covering ap- 
proximately 50% of the stellar surfa ce (if the spots ar e 
cooler than the photosphere by 15%: lRibas et al.ll2008D . 
GJ 3236 shows out of eclipse modulations of comparable 
amplitude to YY Gem, a compatible radius discrepancy 
within the present observational errors, and similarly has 
effective temperatures somewhat cooler than the models 
predict, so this is an attractive hypothesis. By obtain- 
ing precise, multi-band photometry covering the entire 
orbital phase, it may be possible to constrain the spot 
temperatures, and hence observationally test this argu- 
ment. 

The main limitations in the orbital and geometric mod- 
eling of the system arise from the lack of out-of-eclipse 
data in multiple passbands, which leads to large errors 
(4 — 5%) on the radii when we take the degeneracies in 
the spot configuration into account. The precision of 
the mass estimates (5%) is limited by the error in the 
radial velocity measurements, and an important contri- 
bution to the uncertainties in the radii will be from the 
spectroscopic light ratio if the uncertainty in spot pa- 



rameters can be resolved. Further high-precision photo- 
metric measurements and high-resolution spectroscopy 
are therefore needed, and combined with careful analysis 
these should allow the precision of the mass and radius 
measurements for this system to be improved, potentially 
to beyond the 2% level as for the well-known systems 
CM Dra and YY Gem. 
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